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■ Abstract 

We quantize a flat FRW cosmology in the context of the f{R) gravity by Noether symme- 

■ try approach. We explicitly calculate the form of f{R) for which such symmetries exist. It is 
shown that the existence of a Noether symmetry yields a general solution of the Wheclcr-DeWitt 

O , . equation where can be expressed as a superposition of states of the form e''^ . In terms of Hartle 

criterion, this type of wave function exhibits classical correlations, i.e. the emergent of classical 



, universe is expected due to the oscillating behavior of the solutions of Wheeler-DeWitt equation. 



oo 

o 



According to this interpretation we also provide the Noether symmetric classical solutions of our 
f{R) cosmological model. 



PACS numbers: 98.80.Qc, 04.50.+h, 04.20.Fy 
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tJ- ; 1 Introduction 

CP I In recent years, modified theories of gravity constructed by adding correction terms to the usual 

Einstein-Hilbert action, have opened a new window to study the accelerated expansion of the universe. 
It has been shown that such correction terms could give rise to accelerating solutions of the field 
equations without having to invoke concepts such as dark energy [1]. In a more general setting, one 
^ ■ can use a generic function f{R), instead of the usual Ricci scalar R as the action of general relativity. 

Such f{R) gravity theories have been extensively studied in the literature over the past few years, see 
[2] for a review. In finding the dynamical equations of motion one can vary the action with respect to 
the metric (metric formalism) , or view the metric and connections as independent dynamical variables 
and vary the action with respect to both independently (Palatini formalism) [3]. In this theory, the 
Palatini form of the action is shown to be equivalent to a scalar-tensor type theory from which the 
scalar field kinetic energy is absent. This is achieved by introducing a conformal transformation in 
which the conformal factor is taken as an auxiliary scalar field [4]. As is well known, in the usual 
Einstein-Hilbert action these two approaches give the same field equations. However, in f{R) gravity 
the Palatini formalism leads to different dynamical equations due to nonlinear terms in the action. 
There is also a third version of f{R) gravity in which the Lagrangian of the matter depends on the 
connections of the metric (metric-affine formalism) [5]. 

In a previous work [6] , we studied a flat FRW space-time in the framework of the metric formalism 
of f{R) gravity. In [6] we constructed an effective Lagrangian in the minisuperspace {a,R} where a 
and R being the scale factor and Ricci scalar respectively. The form of the function f{R) appearing 
in the modified action is then found by demanding that the Lagrangian admits the desired Noether 
symmetry [7]. A similar study of this issue in the Palatini framework can be found in [8]. By the 



*email: b-vakili@sbu.ac.ir, bvakili45@gmail.com 



1 



Noether symmetry of a given minisuperspace cosmological model we mean that there exists a vector 
field X, as the infinitesimal generator of the symmetry on the tangent space of the configuration 
space such that the Lie derivative of the Lagrangian with respect to this vector field vanishes. For 
some applications of the Noether symmetry approach in various cosmological models see [9] . 

As mentioned above, although the corrections to the results of standard general relativity are 
widely investigated in literature in the f{R) gravity context, these works are often in the classical 
regimes [10]. The cases dealing with quantum f{R) models have seldom been studied in the literature 
[11], and it would be of interest to employ such models in this study. 

In this paper we consider the same model as in [6] and try to quantize it by Noether symmetry 
approach. In general, the existence of a symmetry results in appearing of constants of motion which 
are related to the existence of cyclic variables in the dynamics [12]. This is the key result in our 
quantization procedure. Indeed, in terms of the cyclic variables the conserved quantities are nothing 
but the corresponding conjugate momenta. We shall see that applying the quantum version of these 
symmetries on the wave function of the universe (which satisfies the Wheeler-DeWitt equation) yields 
an oscillatory behavior for the wave function in the direction of the symmetries. On the other hand, in 
the scmiclassical approximation for the canonical quantum gravity, one can show that a wave function 
with classical correlations is a superposition of states of the form e*'^, i.e. with oscillatory behavior 
[13]. With this interpretation scheme for the solutions of Wheeler-DeWitt equation, we also obtain 
the classical trajectory and show that the cosmological scale factor obeys a power law expansion. 

It is to be noted that our presentation of quantum cosmology is applying quantum mechanics 
to a reduced dynamical system, the so-called minisuperspace. This is because that the fundamental 
equation of quantum cosmology, the Wheeler-DeWitt equation, is a differential equation on the 
infinite dimensional superspace and dealing with its solutions in such a space is not an easy task. 
Therefore, in order to solve this equation and obtain the wave function of the universe we can use an 
approximation method in which one truncates the infinite degrees of freedom to a finite dimensional 
submanifold called minisuperspace. In an alternative approach to deal with quantum effects in the 
early universe one can use the idea of quantum field theory in curved space time [14]. In this theory 
we assume that the matter fields are quantized while the gravitational field is classical and given by 
Einstein or any other modified gravity field equations. In this scheme quantum gravity is truncated 
at one-loop level and for computation of the one-loop corrections quantum fluctuations are imposed 
on a given classical background. Our aim in this letter is to study of quantization procedure in a 
Noether symmetric f{R) gravity model in a reduced phase space framework, i.e. as mentioned above 
we first fix the background and then try to quantize it by Wheeler-DeWitt approach. In [15] the role 
of quantum effects to the background cosmology at one-loop level are investigated in f{R) gravity 
framework. 

2 The Noether symmetric phase space of the model 

In this section we consider a spatially flat FRW cosmology within the framework of f{R) gravity. 
Since our goal is to study models which exhibit Noether symmetry, we do not include any matter 
contribution in the action. Let us start from the action (we work in units where c = h = IBttG = 1) 



where R is the scalar curvature and f{R) is an arbitrary function of R. By varying the above action 
with respect to metric we obtain the equation of motion as 



where a prime represents differentiation with respect to R. We assume that the geometry of space- 
time is described by the flat FRW metric which seems to be consistent with the present cosmological 





-5^,/(i?) - i?^./'(i?) + V^V./'(i?) - 9^,,af{R) = 0, 
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observations 

3 
i=l 

With this background geometry the field equations read 



2H + 3H^ = -j 



f"'R^ + 2HRf" + f'R +^if- Rf) 



(4) 



{f'R-f)-6RHf"], (5) 



where H = a /a is the Hubble parameter and a dot represents differentiation with respect to t. 
To study the symmetries of the minisuperspace under consideration, we need an effective point- 
like Lagrangian for the model whose variation with respect to its dynamical variables yields the 
correct equations of motion. Following [16], we consider the action described above as representing 
a dynamical system in which the scale factor a and scalar curvature R play the role of independent 
dynamical variables. In [6, 12, 16] it is shown that such point-like Lagrangian takes the form 

£(a, d, R, R) = Qd^af + QdRa^f" + a^if'R - /). (6) 

The Hamiltonian corresponding to Lagrangian (6) can then be written in terms of a, d, R and R as 

n{a, d, R, R) = Gd'^af + GdRa^f" - a^{f'R - /). (7) 

Therefore, our cosmological setting is equivalent to a dynamical system where the phase space is 
spanned by {a, R,Pa,PR} with Lagrangian (6) describing the dynamics with respect to time t. Now, 
it is easy to see that variation of Lagrangian (6) with respect to R gives the well-known relation for 
the scalar curvature, while variation with respect to a yields the field equation (4). Also, equation 
(5) is nothing but the zero energy condition TC = (Hamiltonian constraint). 

As is well known, Noether symmetry approach is a powerful tool in finding the solution to a 
given Lagrangian, including the one presented above. In this approach, one is concerned with finding 
the cyclic variables related to conserved quantities and consequently reducing the dynamics of the 
system to a manageable one. The investigation of Noether symmetry in the model presented above 
is therefore the goal we shall pursue here. Following [7], we define the Noether symmetry induced 
on the model by a vector field X on the tangent space TQ = (a, R, d, R^ of the configuration space 
Q = (a, R) of Lagrangian (6) as 

d d da d dp d 

such that the Lie derivative of the Lagrangian with respect to this vector field vanishes 

LxC = 0. (9) 

In (8), a and /? are some unknown functions of a and R. Now, it is easy to see that the constants of 
motion corresponding to such a symmetry are [6, 12] 

Q = apa + /3pR. (10) 

In order to obtain the functions a and /? we use equation (9). In general this equation gives a 
quadratic polynomial in terms of d and R with coefficients being partial derivatives of a and P with 
respect to the configuration variables a and R. Thus, the resulting expression is identically equal to 
zero if and only if these coefficients are zero. This leads to a system of partial differential equations 
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for a and /3 [6]- [9]. For the model at hand, this system of differential equations has been investigated 
carefully in [6] and it has been shown that in the case where /" 7^ 0, we have the following solutions 

a{a) = /3ia,R) = Ra''^, (11) 

and ^ 

f{R) = Ri. (12) 



Substituting this results into equation (6) we obtain the Lagrangian of the Noether symmetric model 

(13) 



as 

C = 9a^aR'/^ + ^aRa^R-'/^ + ^a^R^/^ 



The momenta conjugate to variables a and R are 

d£ 9 



p^= = iSaaR^/^ + -a^RR-^/^, (14) 
da 2 

pj, = ^ = %''R-y^ (15) 

dR 2 ^ ^ 

Therefore, the corresponding Hamiltonian reads 

n = la-'R'/'paPR - la-^R'^Vn - ^a^R'/'- (16) 

Although, the classical equations of motion resulting from the Lagrangian (13) or Hamiltonian (16) 
can be solved to give the corresponding classical cosmology, Hamiltonian (16) has not the desired 
form for the construction of the Wheeler-DeWitt equation describing the relevant quantum cosmology. 
Furthermore, the Lagrangian (13) does not exhibit the existence of a cyclic variable corresponded 
to the Noether symmetry. To be more precise, we seek a point transformation (o, R) — > (u, v) on 
the vector field (8) such that in terms of the new variables {u, v), the Lagrangian includes one cyclic 
variable. A general discussion of this issue can be found in [12]. Under such point transformation it 
is easy to show that the vector field (8) takes the form 

X = + + + (17) 

One can show that if X is a Noether symmetry of the Lagrangian, X has also this property, that is 

XC = 0^ XC = 0. (18) 

Thus, if we demand 

Xu = 1, Xv = 0, (19) 

we get 



remark about this form of f{R) function is that it is found only by demanding that the Lagrangian admits the 
Noether symmetry. Although, wc shall see that this form of f{R) gravity yields a power law inflationary cosmology, 
(see [17] for more general f{R) models which lead to the inflationary and late-time accelerating epochs) but in view 
of having the correct wcak-fleld limit at Newtonian and post-Newtonian levels has not a desired form. The conditions 
under which a modified gravity model passes the local and astrophysical tests such as Newton law and solar system 
tests are investigated in [18]. In these works such f{R) theories are studied which satisfy the conditions 

lim f{R) = cons., lim f{R) = 0, 

and shown that they pass Newton law, stability of Earth-like gravitational solution, heavy mass for additional scalar 
degree of freedom, etc.. Equation (12) shows that our Noether symmetric model does not satisfy the above conditions 
and hence is not a viable theory with correct Newtonian and post-Newtonian limits. This is not surprising since it is 
well-known that a large class of f{R) theories suffer from this issue [19]. 
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This means that ii is a cychc variable and the dynamics can be reduced. On the other hand, the 
constant of motion Q becomes 

_ ^ dC ^dC 

Q = apa + Ppr = + fi—^ = 

da dR 

fdCdu dC dv dCdii dCdv\ / dC du dC dv dC dix dC dv\ 
Since (a, R) — >■ {u, v) is a point transformation, we have 



du dv du dv ^ 
da~ da~ dR~ dR~ ' 



dii du du du dv dv dii dv 
dd^da' 'dR^dR' dd^da' 'dR^dR' 



dC du dC dv\ / dC du dC dv \ 



and 



Therefore, 

^ ^dii da^ dv da) ^ ^ \ du dR ' dv dR J 
du ^ ^du\ dC ^ / dv ^ ^dv\dC 
da dRJ dii \ da dRj dv 

.dC ,dC dC , , 

Thus, as expected the constant of motion which corresponds to the Noether symmetry is nothing but 
the momentum conjugated to the cyclic variable. To find the explicit form of the above mentioned 
point transformation we should solve the equations (19), which give 

1 -idu ^ r, du , , 

^ ^dv ^ dv ^ , . , 

These differential equations admit the following general solutions 

u{a, R) = -a^ + (?!)i(ai?^/2), v{a, R) = (j)2{aR^/'^), (25) 

where 0i and 02 are two arbitrary functions of aR^^'^. As is indicated in [12], "the change of coordi- 
nates is not unique and a clever choice is always important." With a glance at the Lagrangian (13), 
we choose the functions (pi and 02 as 

0i(ai?V2) ^ (ai?V2)M^ 02(ai?^/2) ^ ^aR^/^y, (26) 
where /j, and u are some constants. With this choice, the Lagrangian (13) takes the form 

--hrv — ^ — v^ v~uv + -V^'". (27) 

2z^2 2i/ 2 ^ ^ 



It is clear from this Lagrangian that u is cyclic and the Noether symmetry is given hy Pu = Q = cons. . 
Also, the momenta conjugate to u and v are 

dC 9 1^ . dC ^ a ^-21^+1 . 9 hziL . 

Pu = ^ = -TT'" " ^' Pv = ^ = Q^-y - V - —V - u, (28) 
du 2u ov 2u 

which give rise to the following Hamiltonian for our dynamical system 

10 M~lo 2 ly — l 1q/ 

n = — Q-l^v~p^ - -vv~PuPv - 2^ • (29) 

The preliminary set-up for writing the action is now complete. In the next section, we shall focus 
attention on the study of the quantum cosmology of the model described above. 
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3 Quantization of the model 



Standard cosmological models based on classical general relativity have no convincing precise answer 
for the presence of the so-called Big-Bang singularity. Any hope of dealing with such singularities 
would be in vein unless a reliable quantum theory of gravity can be constructed. In the absence of 
a full theory of quantum gravity, it would be useful to describe the quantum states of the universe 
within the context of quantum cosmology [20]. In this formalism which is based on the canonical 
quantization procedure, one first freezes a large number of degrees of freedom and then quantizes 
the remaining ones. The quantum state of the universe is then described by a wave function in the 
minisuperspace, a function of the 3-geometry of the model and matter fields presented in the theory, 
satisfying the Wheeler-DeWitt equation, that is, Ti.^ = 0, where H is the operator form of the 
Hamiltonian given by equation (29) and ^ is the wave function of the universe. On the other hand, 
the existence of a Noether symmetry in the model reduces the dynamics through Pu = Q, where its 
quantum version can be considered as a constraint = Q^. Therefore, the quantum cosmology 
of our Noether symmetric f{R) model can be described by the following equations 



— /ii; - Pu - PvV +v pyV )Pu- -^v ' 



'i!{u,v) = 0, 



(30) 



Pu^{u,v) = Q^{u,v), (31) 

where the parameters r and s satisfy r + s = and denote the ambiguity in the ordering of factors 
V and Pv in the second term of (29). With the replacement Pu — — and similarly for p^ the above 
equations read 
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+ 



1 



du 3 dudv 2 



d 

-i—^{u,v) = Q^{u,v). 



'^{u,v) = 0, 



(32) 



(33) 



The solutions of the above differential equations are separable and may be written in the form 
^{u,v) = U{u)V{v). Equation (33) can be immediately integrated leading to a oscillatory behavior 
for the wave function in u direction, i.e. in the direction of symmetry, that is 



iQu 



Viv). 



(34) 



Substitution this result into relation (32) yields the following equation for the function V{v) 

3 . 4^ v — 1 



dV 
dv 



-—i—Uv V 
3 AuQ 



IV V 



2v 



V{v), 



with solution 



V [v) = v 2i/ e ^^'3 '. 

Thus, the eigenfunctions of the equations (32) and (33) can be written as 

^nu{u,v) = V e ^ 3^ 16Q 



(35) 
(36) 

(37) 



We may now write the general solutions to the Wheeler-DeWitt equation as a superposition of the 
eigenfunctions with shifted Gaussian weight functions 



J —oo J—oo 



(38) 



We see that the wave function is a superposition of states of the form e^^ . In scmiclassical ap- 
proximation for quantum gravity [13], this type of state represents the correlations between classical 
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trajectories and the peaks of the wave function [12]. Inserting ^ e**^ into Wheeler-DeWitt equation, 
we are led to the Hamilton-Jacobi equation for S. Thus, the classical trajectories can be obtained 
by rewriting the momenta as derivative of S with respect to the corresponding variables, that is, 
Pq = Therefore, in semiclassical limit, by identifying the exponential factor of (37) with S, we 
can recover the corresponding classical cosmology 

The classical trajectories, which determine the behavior of the scale factor and Ricci scalar are given 

by 

_dS_ _dS_ 

Pu — ) Pv — 

OU OV 

Using the definition of Pu and py in (28), the equations for the classical trajectories become 

-2;;^~^ = (40) 

a H-2U+1 9 inii . 5 /X 3 4^ 

" V V " u = —-—Qv ^ — " • (41) 

Equation (40) can be easily integrated leading to 

v{t) = (Qt - tor , (42) 

where Q = — 1(5 aiid to is an integrating constant. Substituting the above results into equation (41) 
yields 



u{t) = l{Qt- tor - (Qt - to)' + C, (43) 



where C is another integrating constant which we can choose it to be zero. Going back to the variables 
a and R, we obtain the corresponding classical cosmology as 



2 



«(*)= :^{Qt-tor + Y^iQt-toY 

Qt - toY 



1/2 



(44) 



R(t) = T. (45) 

liQt - tor + Toh (Qt - to)' 

For t to/Q, we have an initial singularity in which the scale factor goes to zero while the Ricci 
scalar has a large value. On the other hand, in the late time, the universe evolves with a power law 
expansion {a{t) ~ for < /x < 4 and a{t) ~ t'*/^ for /x > 4) and the scalar curvature goes to zero 
in this limit. 

In general, one of the most important features in quantum cosmology is the recovery of classical 
cosmology from the corresponding quantum model, or in other words, how can the Wheeler-DeWitt 
wave functions predict a classical universe. We see that the oscillatory solutions of the form * ~ e*'^ 
for Wheeler-DeWitt equation yield the classical solution (44) which for positive values of n can be 
viewed as an accelerating cosmology. Since the quantum effects in cosmology are important in the 
very early times of cosmic evolution, we immediately see that, in this limit, the scale factor has the 
behavior a{t) ~ t^^^ for < < 4 and a{t) ~ for /x > 4, which for /x > 2 describe a power law 
inflationary behavior. The universe then undergoes to a late-time accelerating phase, also with a 
power law expansion behavior. An important question in inflationary models is how much inflation 
do the model predict or in other words, what is the mechanism through which the universe exits from 
the inflationary epoch and undergoes into radiation or matter-dominated eras. As is well known, this 
is largely depend on the behavior of a scalar field with which the universe nucleates [21]. In f{R) 
gravity models some of such mechanisms describing the transition between different epochs of cosmic 
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20.5 



Figure 1: The square of wave function. The figure is plotted for numerical values ? = 10, a = 5, = 3, z-'o = and 



evolution are proposed in [22] by introducing a scalar field cf). Indeed, it is easy to verify that the 
metric f{R) gravity action (1) is dynamically equivalent to [2] 



where (/> = f'{R) and ¥{4)) = (f>R{^) - f{R{(t>))- In our f{R) ~ i?3/2 model we have ~ (f)'^ and 
hence V{4)) ~ 4>^ . The dynamics of this scalar field may show the the possibility of the transition 
between different eras of cosmic evolution. But as is indicated in [22] "we should note that in this 
case, there is no matter and /(-R)-terms contribution plays the role of the matter instead of the real 
matter". Therefore, it is to be noted that our presentation does not claim to clear the role of inflation 
scenario in a fundamental way because we just study the problem in a special simple model. However, 
this may reflect realistic scenarios in similar investigations which deal with this problem in a more 
fundamental way [22]. 

From equations (42) and (43) we see that the classical trajectories obey the relation 



As we have mentioned above (see (38)), we are looking for a coherent wave function with a good 
asymptotic behavior in the minisuperspace and peaking in the vicinity of the classical loci (47) in the 
configuration space spanned by {u,v}. It is well known that the general solution of Wheeler-DeWitt 
equation may be constructed by superposition of its eigenfunctions which in our problem at hand 
are labeled by /i and v. Therefore, the wavepacket (38) is what we need. We take the solution as 
being represented by equation (38) with the integrals to be truncated at suitable values of // and u 
displaying this peak. Figure 1 shows the square of wave function for typical values of the parameters 
where we have taken the integrals from to 2 for // and from -15 to 15 for v. It is seen that almost 
a good correlation exists between this pattern and the classical trajectories. 

4 Conclusions 

In this letter we have studied a generic f{R) cosmological model by Noether symmetry approach. For 
the background geometry, we have considered a flat FRW metric and derived the general equations 
of motion in this background. The phase space was then constructed by taking the scale factor a 
and Ricci scalar R as the independent dynamical variables. The Lagrangian of the model in the 



Q = 3. 




(46) 



configuration space spanned by {a, R} is so constructed such that its variation with respect to these 
dynamical variables yields the correct field equations. The existence of Noether symmetry implies 
that the Lie derivative of this Lagrangian with respect to the infinitesimal generator of the desired 
symmetry vanishes. In [6] we have shown that by applying this condition to the Lagrangian of the 
model, one can obtain the explicit form of the corresponding f{R) function which has led us to the 
Lagrangian (13) and Hamiltonian (16) of the Noether symmetric model. Since the Lagrangian (13) 
does not exhibit the existence of a cyclic variable corresponded to the Noether symmetry, we have 
provided a point transformation {a,R) {u,v) such that in terms of the new variables {u,v), the 
Lagrangian includes one cyclic variable. We have then quantized the model and shown that the 
corresponding quantum cosmology and the ensuing Wheeler-DeWitt equation are amenable to exact 
solutions in terms of a superposition of states of the form e*'^ due to the existence of Noether symmetry. 
In semiclassical approximation for quantum gravity, this type of state represents the correlations 
between classical trajectories and the peaks of the wave function. Using this interpretation we have 
shown that the corresponding classical cosmology results in a power law accelerated expansion for 
the scale factor of the universe either in its early or late time evolution. 
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